Calculation of P, T-odd effects in 205 T1F including electron correlation 
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A method and codes for two-step correlation calculation of heavy-atom molecules have been 
developed, employing the generalized relativistic effective core potential (GRECP) and relativistic 
coupled cluster (RCC) methods at the first step, followed by nonvariational one-center restoration 
of proper four-component spinors in the heavy cores. Electron correlation is included for the first 
time in an ab initio calculation of the interaction of the permanent P, T-odd proton electric dipole 
moment with the internal electromagnetic field in a molecule. Inclusion of electron correlation by 
GRECP/RCC has a major effect on the P, T-odd parameters of 205 T1F, decreasing M by 17% and 
X by 22%. 

PACS numbers: 31.90.+S, 35.10.Wb, 31.20.-d 



Introduction. The measurement of permanent electric 
dipole moments (EDM) of elementary particles is highly 
important for the theory of P, T-odd interactions. Exper- 
iments performed so far have given only upper bounds for 
the EDMs. The extraction of EDMs from measurements 
on molecules containing heavy atoms requires knowledge 
of nuclear and electronic properties of the molecule. High 
quality calculations of the relevant electronic properties 
are therefore essential for accurate determination of the 
EDMs 0| . 

Here we consider the interaction of the proton EDM 
with the internal electromagnetic field of the 205 TIF 
molecule. This molecule is one of the best candidates for 
proton EDM measurements. Following Hinds and San- 
dars m, the effective interaction with the proton EDM 
in TIF is written in the form 



H cB = {d v + d M )a N • A , 



(1) 



where $n is the Tl nuclear spin operator, A is the unit 
vector along the internuclear axis z from Tl to F, d v and 
d M are constants corresponding to the volume and mag- 
netic effects according to Setoff's theory 0. Hinds and 
Sandars showed || that the volume effect in a coordinate 
system centered on the Tl nucleus is given by 
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where d p is the proton EDM, R is a factor determined 
by the nuclear structure of 205 T1, 
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Pip(r) is the electronic density calculated from the wave 
function ip. Keeping only the dominant diagonal terms 
of the two-electron operator for the magnetic effect (see 
U) they have 



d M = lSd v 



2mc 



M , 



(4) 



where /i, m and Z are the magnetic moment, mass and 
charge of the Tl nucleus, c is the velocity of light, 
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li is the orbital momentum operator of electron i, and cti 
are its Dirac matrices. H c g leads to different hyperfine 
splitting of TIF in parallel and antiparallel electric and 
magnetic fields. The level shift hv = 4(d v +d M )(aN-^) is 
measured experimentally (for the latest data see f|; an- 
other experiment is now in preparation at the Petersburg 
Nuclear Physics Institute). 

The parameters X of Eq. (||) and M of Eq. (||) are 
determined by the electronic structure of the molecule. 
They were calculated recently for the X0 + ground state 
of TIF by Parpia || and by Quiney et al. using the 
Dirac-Hartree-Fock (DHF) method with large Gaussian 
basis sets (see Table ||). No calculation which includes 
correlation effects is available. The main goal of the 
present work is to calculate the X and M parameters 
for the molecule with correlation included to high order. 

Methods. The generalized relativistic effective core 
potential (GRECP) method § is applied to the TIF 
molecule. A two-component electronic (pseudo)wave 
function is first obtained with the 21-electron GRECP 
H [l^] for Tl, providing proper electronic density in the 
valence and outer core regions, followed by restoration 
of the proper shape of the four-component molecular 
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spinors in the inner core region of Tl. Details of the 
method may be found elsewhere [|b[ [Tl[ [l2[ [l3|. 

The correlation spin-orbital basis set used consisted of 
26s, 25p, 18c?, 12/, and lOg Gaussian-type orbitals on Tl, 
contracted to 6s6p4:d2flg. The basis was optimized in 
a series of atomic two-component GRECP calculations, 
with correlation included by the all-order relativistic cou- 
pled cluster (RCC) method Q with single and double 
excitations; the average energy of the two lowest states of 
the atom was minimized. The basis set generation pro- 
cedure is described in Refs. |l5|, |l(|. The basis set was 
designed to describe correlation in the outer core 5s and 
5p shells of Tl, in addition to the 5d and valence shells. 
While 5s and bp correlation may not be important for 
many of the chemical and physical properties of the atom, 
it is essential for describing properties coming from inner 
regions, including P, T-odd effects. The (14s9p4d3/)/- 
[<is3p2dlf] basis set from the ANO-L library |l7| is used 
for fluorine. 

A one-component self-consistent-field (SCF) calcula- 
tion of the (lcr. . .7cr) 14 (l7r27r37r) 12 (l<5) 4 ground state 
of TIF is performed first, using the GRECP for Tl 
which simulates the interactions of the valence and 
outer core (bsbpbd) electrons with the inner core 
[Kr]4d 4 /2 4d6 /2 4/f /2 4/8 /2 . This is followed by two- 
component RCC calculations, with only single (RCC- 
S) or with single and double (RCC-SD) cluster ampli- 
tudes. The RCC-S calculations with the spin-dependent 
GRECP operator take into account effects of spin-orbit 
interaction at the level of the one-configurational SCF- 
type method. The RCC-SD calculations include, in ad- 
dition, the most important electron correlation effects. 

The electron density obtained from the two-component 
GRECP/RCC (pseudo)wave function in the valence and 
outer core regions is very close to that of the correspond- 
ing all-electron four-component function. The pseu- 
dospinors are smoothed in the inner core region |8J, so 
that the electronic density in this region is not correct. 
The operators in equations (|^) and (^|) are heavily con- 
centrated near the nucleus, and are therefore strongly 
affected by the wave function in the inner region. The 
four-component molecular spinors must therefore be re- 
stored in the inner region of Tl. All molecular spinors 
fa are restored as one-center expansions on the Tl nu- 
cleus, using the nonvariational restoration scheme (see 
||, [ll], |l2|, [13| and references therein). 

The restoration is started by generating equivalent ba- 
sis sets of atomic (one-center) four-component spinors 

fnlj { r )Xljm 

g n ij(r)xi'jm 

{fnij{r)xijm} by atomic finite-difference all-electron 
DHF and two-component GRECP/SCF calculations of 
the same valence configurations of Tl and its ions. Here 
n is the principal quantum number, j and m are the total 
electronic momentum and its projection on the internu- 
clear axis, I and I' are the orbital momenta, and l'=2j—l. 
The nucleus is modeled as a uniform charge distribution 
within a sphere with radius r nuc i = 7.1fm = 1.34 x 10~ 4 



and two-component pseudospinors 



a.u., whereas previous calculations employed a spherical 
Gaussian nuclear charge distribution |], [?J (the root mean 
square radius in all calculations is 5.5 fm, in accord with 
the parametrization of Johnson and Soff [fl8| , and agrees 
with the experimental value 5.483 fm for the 205 T1 nu- 
cleus Jl^] ) . The all-electron four-component hfd and 
two-component GRECP / hfj || [TXJ codes were employed 
to generate the two equivalent [15sl2pl2<f8/] numerical 
basis sets for restoration. These sets, describing mainly 
the core region, are generated independently of the ba- 
sis set for the molecular GRECP calculations discussed 
earlier. The molecular pseudospinorbitals are then ex- 
panded in the basis set of one-center two-component 
atomic pseudospinors, 
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Note that for linear molecules only one value of m sur- 
vives in the sum for every fa. Finally, the two-component 
pseudospinors in the basis are replaced by the equivalent 
four-component spinors and the expansion coefficients 
from Eq. (§) are preserved [□} [TJ, [TJ]: 



L^] = \l + l/2\ 

f- ^ I ^ nljm \9nij(r)x 
1=0 j=\l-l/2\ n,m v 
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The molecular four-component spinors constructed this 
way are orthogonal to the inner core spinors of Tl, as the 
atomic basis functions used in Eq. (0) are generated with 
the inner core electrons treated as frozen. 

The quality of the approximation for the two-center 
molecular spinors and, consequently, of the calculated 
properties increases with the value of I/ max . A series of 
calculations of the M parameter was performed using Eq. 
(0) with basis functions going up to p, d and / harmonics. 
We found (see Table |) that including only s and p func- 
tions in the expansion determines M with 90% accuracy. 
Because the contribution of / is only about 0.3% and 
amplitudes of higher harmonics on the nucleus are sup- 
pressed by the leading term ~ pO' - V 2 ^ the error due to 
the neglect of spherical harmonics beyond / is estimated 
to be below 0.1%. Calculation of the X parameter re- 
quires s and p harmonics (see Ref. 0), although, strictly 
speaking, d harmonics also give nonzero contributions. 

The restoration procedure implemented here gives a 
very good description of the wave function in the core 
region, which is important for accurate evaluation of the 
X and M parameters. This is done at a fraction of the 
cost necessary for all-electron four-component molecu- 
lar calculations with Gaussian basis sets, where a large 
number of additional basis functions must be included 
for proper description of the inner core region and small 
components of spinors j7|. Here we calculate (restore) 
the four-component electronic wave function in the core 
region from the (pseudo)wave function obtained in the 
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molecular GRECP calculation, which may be considered 
"frozen" in the valence region at the restoration stage. 
Basis functions describing a large number of chemically 
inert core electrons may thus be excluded from the molec- 
ular GRECP calculation. 

The X and M parameters were calculated by the fi- 
nite field method (see, e.g., Refs. |2l], |2^|). The operator 
corresponding to a desired property [see Eqs. (||) and (||)] 
is multiplied by a small parameter A and added to the 
Hamiltonian. The derivative of the energy with respect 
to A gives the computed property. This is strictly correct 
only at the limit of vanishing A, but it is usually possible 
to find a range of A values where the energy is linear with 
A and energy changes are large enough to allow sufficient 
precision. The quadratic dependence on A is eliminated 
in the present calculations by averaging absolute energy 
changes obtained with A of opposite signs. 

Results and discussion. Calculations were carried out 
at two internuclear separations, the equilibrium R e = 
2.0844 A as in Ref. ||, and 2.1 A, for comparison with 
Ref. @. The results are collected in Table | The first 
point to notice is the difference between spin-averaged 
SCF and RCC-S values, which include spin-orbit interac- 
tion effects. These effects increase X by 9% and decrease 
M by 21%. The RCC-S function may be written as a sin- 
gle determinant, and results may therefore be compared 
with DHF values, even though the RCC-S function is not 
variational. GRECP /RCC-S values of the M parameter 
are indeed within 3% and 1% of corresponding DHF val- 
ues H (Table ||). This agreement confirms the validity 
of the approximations made by us. In particular, freezing 
the inner core shells is justified, as inner core relaxation 
effects have little influence on the properties calculated 
here, a conclusion already drawn by Quiney et al. @. 

Much larger differences occur for the X parameter. 
Here there are also large differences between the two DHF 
calculations, which cannot be explained by the small 
change in internuclear separation. The value of X may 
be expected to be less stable than M, because it is deter- 
mined by the derivative of the electronic density at the Tl 
nucleus and involves large cancellations Q between con- 
tributions of large and small components, each of which 
is about 20 times larger than their sum. Thus, a strong 
dependence of X on the basis used may be expected. 
The DHF values collected in Table | indeed show such 
dependence. Results obtained in Refs. || and M with 
comparable even-tempered basis sets, (28s28pl2c?8/) and 
(28s28pl4d8/), are rather close, differing by 340 a.u. Im- 
proving the Tl basis to (34s34pl6d9/) M increases X by 
650 a.u. or 8%. Further improvement of the basis may 
be expected to yield even higher X values. The numeri- 
cal basis functions obtained in atomic DHF calculations 
and used for the restoration are highly accurate near the 
nucleus, so that our RCC-S value for X, which is higher 
than that of Quiney et al. 0, seems reasonable. The dif- 
ferent nuclear models used in the present and DHF [|[ 
calculations may also contribute to the disagreement in 



X, which is determined by the derivative of the electronic 
charge density at a single point, the Tl origin. M is af- 
fected by ip in a broader region, and is therefore far less 
sensitive to the nuclear model. 

The main goal of this work is the evaluation of electron 
correlation effects on the P, T-odd parameters. These ef- 
fects are calculated by the RCC-SD method at the molec- 
ular equilibrium separation R e . A major correlation con- 
tribution is observed, decreasing M by 17% and X by 
22%. 

Using the correlated values for X and M calculated 
here and R = 1.036 x 10 -9 a.u. from ^3|, one obtains 
from Eqs. @ and ([§) 

d v = -7.909 x 10"% a.u. (8) 

d M = 1.622 x 10~% a.u. (9) 

The effective electric field interacting with the EDM of 
the valence proton of 205 Tl in the fully polarized TIF 
molecule is E = \d v + d M \/d p = 6.287xl0~ 6 a.u. = 32.33 
kV/cm; the revised proton EDM limit for the experiment 
of Ref. § is d p = (-1.7 ± 2.8)xl0~ 23 e-cm. 

The hyperfine structure constants of Tl Qp\/ 2 an d Tl 2+ 
6s 1 , which (like X and M) depend on operators concen- 
trated near the Tl nucleus, were also calculated. The 
errors in the DF values are 10-15%; RCC-SD results are 
within 1-4% of experiment. The improvement in X and 
M upon inclusion of correlation is expected to be similar. 

Concluding remarks. Note that the codes developed 
for GRECP/RCC calculation followed by nonvariational 
one-center restoration in heavy cores are equally appli- 
cable to calculation of other properties described by op- 
erators singular near nuclei (hyperfine structure, quan- 
tum electodynamic effects, etc.). Because the Fock-space 
RCC-SD approach (l4) is used, the two-step method is 
applicable to both closed-shell and open-shell systems, 
including excited states. In particular, calculations for 
the ground state of YbF and for excited states of PbO 
are in progress now. Triple and higher cluster ampli- 
tudes in the valence region are important for chemical 
and spectroscopic properties, but not for the effects dis- 
cussed here, as concluded from previous calculations for 
YbF fl2]| . These excitations are believed to be unim- 
portant in the core region too. We therefore suggest that 
further improvement in the correlation treatment will not 
seriously affect our M and X values. 
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TABLE I: Calculated X and M parameters [Eqs. (||) and (|)] for the 205 T1F ground state, compared with DHF values with 
different basis sets |], Q . Individual shell contributions are calculated from spin-averaged GRECP /SCF orbitals. GRECP /RCC- 
S results include spin-orbit interaction, and GRECP/RCC-SD values also account for electron correlation. All values in a.u. 







R e 


= 2.0844 A 




R = 2.1 A 


Expansion 


s,p 


s,p,d 


s,p,d,f 


s, p 


s,p 


s,p,d,f 


s, p 


Shell : main contribution 


M 


X 


M 


X 


la 2 : U*{F) 


0.01 


0.02 


0.02 


3 


0.00 


0.02 


1 


2a 2 : 5s 2 (Tl) 


-2.49 


-2.49 


-2.49 


-1114 


-2.44 


-2.44 


-1089 


3cr 2 : 5p 2 (TZ) 


4.21 


3.91 


3.91 


1897 


4.10 


3.82 


1851 


4cr 2 : 2s 2 (F) 


-0.79 


-0.64 


-0.64 


-358 


-0.74 


-0.60 


-335 


5cr 2 : 5d 2 2 (Tl) 


-0.01 


-0.04 


-0.05 


-2 


-0.01 


-0.05 


-2 


6a 2 ■ (6s(T/) + 2p z (F)) 2 


-9.38 


-10.05 


-10.06 


-4414 


-9.38 


-10.02 


-4422 


la 2 : (6s(Tl) - 2p z (F)) 2 


28.13 


27.19 


27.19 


12954 


27.98 


27.07 


12893 


Ik 4 : 5pl5p 2 y (Tl) 


0.00 


-0.26 


-0.26 





0.00 


-0.25 





27T 4 : bdl z bd 2 yz (Tl) 


0.00 


0.31 


0.30 





0.00 


0.27 





37T 4 : 2p 2 2p 2 (F) 


0.00 


-0.39 


-0.40 





0.00 


-0.38 





= 5dl 2 _ v2 Bdl y (Tl) 


0.00 


0.00 


-0.02 





0.00 


-0.02 





Total SCF(spin-avcragcd) 


19.67 


17.56 


17.51 


8967 


19.52 


17.43 


8897 


GRECP/RCC-S 


16.12 




13.84 


9813 


16.02 


13.82 


9726 


DHF ffl Tl:(28s28pl2d8/) 


15.61 


7743 






DHF jpf Tl:(25s25pl2d8/) 
Tl:(28s28pl4d8/) 
Tl:(31s31pl5d8/) 
Tl:(34s34pl6d9/) 






13.64 a 
13.62 a 
13.66 a 
13.63° 


8098 
8089 
8492 
8747 


GRECP/RCC-SD 


11.50 


7635 







a M is calculated in Rcf. pj using two-center molecular spinors, corresponding to infinite L max in Eq. (I 
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